UNCLASSIFIED 


2m6« 


lut  Ute 

ARMED  SERVICES  TECHNICAL  INFORMALON  AGENCY 
ARLLNGTON  HALL  STATION 
ARLINCIIW  12,  VIRGINIA 


UNCLASSIFIED 


MOriCS:  Vhen  fovexnMnt  or  other  drcvloge,  sped* 
flcatloDS  or  other  date  axe  used  for  any  puxpoae 
other  than  In  connection  ^  h  a  definitely  related, 
goremaent  proeureaent  onexetion,  the  U.  S. 
CtoTemnent  thereby  incurs  h-  respoostbility,  nor  any 
obligation  idiatsoever;  and  the  fact  that  th.  C  >'/cxn> 
Bent  SHky  L^ve  foiaulated,  fomiabAd,  cr  xn  any  way 
supplied  the  sa^d  dznwings,  Bpeciflcatians>  or  other 
data  is  not  to  be  regarded  by  iaplication  or  other¬ 
wise  as  in  any  sMnner  licena^sg  the  bolder  or  any 
other  person  or  corporation,  or  conveying  any  rif^ts 
or  pexmission  to  aanufaetuxe,  use  or  sell  any- 
patented  invention  that  may  in  any  way  be  rtlated 
thereto. 


:  i  ■  /■“ 


2256 


APPLICWriOR  Gg  VARUfflONAL  BQDATXOH  OP  MOTIOP 
TO  THE  MONLINBAR  VIBRATXOK  AHALYSXS 
OP 

nO'^OGBNEOl.’S  AID  XAVEnO)  PLATES  AID  SHELLS 


By 


Yl-YU«n  YU 


Air  Porc«  offic*  of  Sciontific  Resuarch 
Contrac*,  AP4r)(63ft)-453 
Technical  Mote  Mo. 14 


POLYTECHNIC  INSTITUTE  OF  BROOKLYN 

DEPARTMENT  OF 
MECHANICAL  ENGINEERING 


Pebruary  1962 


AFOSR  ^  2256 


APPLICATION  OF  VARIATIONAL  EQUATION  OP  >»OTION 
TO  THE  NONLINEAR  VIBRATION  ANALYSIS 
OF 

»«>'^OGENEOrS  AH>  LAYERED  PLATES  AND  SHELLS 


By 

Yi-YUan  Yu 


Air  Force  Office  of  Scientific  Research 
Contract  AP49(63P)-45j 
Technl.?  tl  Note  No.  14 


February  1962 


APPLICATION  OP  VARIATIONAL  ECJl^ATION  OF  mqtION 
TO  THS  NONLINEAR  VIBRATION  ANALYSIS 
OF 

T^O«OCENEOtlS  AND  LAYERED  PLATES  AND  SHELLS^ 
Yi-Yuan  Yu^ 

ABSTRACT 


An  integrated  procedure  is  presented  for  applying 
the  variational  equation  of  motion  to  the  approximate 
analysis  of  nonlinear  vibrations  of  homogeneous  emd 
layered  plates  and  shells  involving  large  deflections. 

The  procedure  consists  of  a  *iequence  of  variational 
approximations.  The  first  of  these  involves  an  approxi¬ 
mation  in  the  thickness  direction  and  yields  i  system  of 
equations  of  motion  and  boundary  conditions  the  plate 
or  shell.  Subsequent  variational  approximati  ‘ns  with  re¬ 
spect  to  the  remaining  space  coordinates  and  ♦•ime,  wher¬ 
ever  needed,  lead  to  a  solution  to  the  nonlir-'ar  ’-ibration 
problem.  The  procedure  is  illustrated  by  a  study  of  the 
nonlinear  free  vibrations  of  homogeneous  and  s  -ndwich  cy¬ 
lindrical  shells,  and  it  appears  to  be  applicable  to  still 
many  other  homogeneous  and  composite  elastic  westerns. 
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INTROPUCTIOM 


In  8oi\'ing  certain  equilibrium  euid  vibration  problems 
in  the  linear  or  nonlinear  theory  of  elasticity,  two  dis¬ 
tinct  and  unrelated  steps  are  often  taken.  In  the  case  o£ 
a  plate  or  shell  problem,  for  instance,  an  approximate 
system  of  equations  that  governs  the  problem  is  usually  de¬ 
rived  first  by  the  use  of  one  of  a  wide  variety  of  avail¬ 
able  methods.  Then,  when  the  system  of  equations  deduced 
cannot  be  solved  exactly^  a  %iride  variety  of  methods  is 
again  available  for  obtaining  an  approximate  solution  of 
the  equations.  The  method  used  for  deriving  the  approxi¬ 
mate  equat:'cns  ar-9  that  used  for  solving  the  equations 
usually  bear  no  relation  to  each  other. 

One  of  the  main  purposes  of  this  paper  is  to  advo¬ 
cate  an  integrated  approximate  procedure  of  solving  a 
large  class  of  problems  in  the  linear  or  nonlinear  theory 
of  elasticity,  and  in  particular,  problems  of  plates  emd 
shells  cf  the  layered  as  %«ell  as  homogeneou^f  type  of  con¬ 
struction.  solely  on  the  basis  of  the  variatX::n«x  equation 
of  motion.  It  is  an  integrated  procedure  in  that  the 
aforementioned  two  steps  are  no  longer  unrelated  to  each 
other.  In  fact,  the  procedure  consists  of  a  sequence  of 
variational  approximations  with  respect  to  the  space  and 
time  coordinates,  carried  out  in  relation  to  the  differ¬ 
ential  equations  and/or  the  boundary  condicicns.  Although 
not  much  originality  can  be  claimed  on  the  variational 
approach,  the  treatment  does  make  the  fullest  systematic 
use  of  the  variational  equation  of  motion.  Besides,  it 
not  only  integrates  some  of  the  variational  approximations 
which  have  been  known  only  as  unrelated  individual  pro¬ 
cedures,  but  it  also  reveals  the  possibility  of  having 
more  general  variational  approximations. 

The  variational  equation  of  motion  in  the  theory  of 
elasticity  is  a  direct  consequence  of  Hamilton's  principle 
and  applies  to  both  linear  and  nonlinear  cases.  The  usual 
formulation  of  the  equation,  as  given  in  Iiove's  book  fll. 
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ooataiM  only  th*  voIwm  and  rarfae*  intagrals  with  re-> 
ayact  to  the  space  coordinates.  In  this  paper  we  shall 
fnrthsr  include  the  integration  with  respect  to  tine  as 
a  necessary  part  in  the  forsnlation  of  the  variational 
equation  of  aotion.  Without  the  additional  tiaw  inte¬ 
gration.  variational  i^proriauitions  can  only  be  perfom- 
ed  with  respect  to  the  apace  coordinates. 

m  idiat  follows  the  proposed  procedure  is  first  out¬ 
lined.  Mien,  as  illustrati«M  and  as  problsaui  of  interest 
by  thesuMlves.  ^proximate  systasHi  of  nonlinear  equations 
of  notiM)  end  boundary  conditions  of  henogensous  and  sand¬ 
wich  cylindrical  shells  are  derived  and  subsequently  solved 
for  the  cases  of  axially  syaewtrical  vibraticne  of  closed 
shells  with  iwovable  hinged  odgsw.  hesults  for  the  non¬ 
linear  frequencies  are  finally  discussed. 
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Tli>  Iat#qr«fd  VMrlatloiial  Approach 

On  tlm  basis  of  tbs  linsar  or  nonlinear  variational 
equation  of  notion  in  a  fom  as  has  just  been  specified, 
the  solutions  to  a  large  class  of  probleau  in  the  theory 
of  elasticity,  and  in  particular  probleu  of  hoeogeneous 
and  layered  plates  and  shells,  nay  be  obtained  by  carry- 
ing  out  a  sequence  of  successive  variational  approxinations. 

In  the  first  of  these  approxinations.  and  in  the  case  of 
plate  and  shell  problens,  for  instance,  the  dependence  of 
the  displacenents  on  the  thidtness  coordinate  is  assumed, 
and  integration  is  carried  out  with  respect  to  this  co¬ 
ordinate.  Th's  firs^  variational  approximation  thus  consti¬ 
tutes  essentially  the  process  of  deriving  plate  or  shell 
equations  or  other  approximate  equations.  Tiowrver,  even 
in  such  a  process  in  this  first  step,  the  varic.tional  equation 
of  motion  does  not  appear  to  have  been  fully  made  use  of 
before,  since,  until  recently,  only  the  volume  integral, 
and  only  that  in  tne  linear  variational  equatim  of  motion 
(without  the  time  integration),  has  been  es^ployed  in  the 
derivation  of  linear  differential  equations  of  motion  of 
plates  [2-7l  and  shells  [C'’.  the  surface  integral  in 
either  the  linear  or  nonlinear  Variational  equation  of  tt 
aiotion  is  believed  to  have  been  used  for  the  first  time 
in  the  recent  derivation  of  the  appropriate  boundary  con¬ 
ditions  in  reference  9  for  sandwich  plates.  Al-hough  the 
boundary  conditions  (as  well  as  the  equations  of  motion) 
may  also  be  derived  by  other  means,  the  integrated  treaty ■ 
ment  on  the  basis  of  the  variation  equation  of  motion  has 
the  advantage  of  being  simple  and  straightforward,  and  it 
permits  the  surface  traction  terms  that  appear  in  the 
boundary  conditions  to  be  incorporated  insmdiately  in  the 
equations  of  motion,  which  is  particularly  desirable  in 
nonlinear  eases. 

Sosmtimes  even  the  approximate  system  of  differential 
equations  and  boundary  conditions  cannot  be  solved  exactly 
for  a  given  equilibrium  nr  vibration  problem.  A  well-known 
approximate  i^oeedure  named  usually  after  Calerkin  may  often 
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be  need,  in  ehieh  dleplneunte  ere  eMuned  such  as  to 
satisfy  ezsetlv  tbs  boundary  conditions  b»t  not  tbs 
differential  equations,  bn  alternative  approxiMte  ihco- 
cedsre  is  that  doe  to  Treffts  [10]  in  idiieb  tbs  assaned 
displacensnts  satisfy  exactly  tbe  differential  equatiems 
but  not  tbe  boundary  conditions.  A  natural  generalisa* 
tion  of  tbe  tMo  appears  to  be  one  in  idiicb  tbe  assusied 
displaceswnts  satisfy  exactly  sane  of  tbe  differential 
equatiems  and  boundary  conditions  and  are  ttiade  to  satis¬ 
fy  tbe  renaining  differential  equations  utd  boundary 
cemditions  approxiamtely  in  tbe  vauriatimtal  sense.  The 
use  of  tbe  generalised  procedure  remains  to  be  explored, 
but,  togetbei  with  its  above  two  variants,  tbe  procedure 
clearly  may  be  considered.  ;ind  in  effect  is.  a  second 
variational  approximation  in  tbe  solution  of  an  elasticity 
problem  baaed  on  tbe  variational  equation  of  motion.  In 
equilibrium  problems  this  is  also  tbe  final  variaticmal 
approximation  that  is  needed. 

For  problems  of  vibration,  and,  in  part1'..-ular,  non¬ 
linear  vibration,  a  third  and  final  variationol  a^roxi- 
mation  with  respect  to  time  is  often  useful.  It  is  this 
last  step  that  needs  tbe  integration  with  respect  to  time 
which  has  been  included  as  part  of  the  variational  e- 
quation  of  motion.  Tbe  approximation  consists  of  essen¬ 
tially  another  ^plication  of  tbe  Galertin  procedure. 

Tbe  successful  use  of  ft  in  solving  nonlineu  /ibration 
problems  involving  single-dsgrse-of-frssdom  systems  has 
been  dsmonstrated  by  Klotter  [ill.  «dx>  prefers  to  call  it 
tbe  hits  Procedure.  In  reference  11  it  is  also  mentioned 
that  tbs  same  procedure  may  be  applied  to  nonlinear  vibra¬ 
tion  analysis  of  two-degrees-of- freedom  systems.  Purtbsr 
applicability  of  tbe  procedure  to  composite  continuous 
systems  has  been  dsmonstrated  in  reference  9  where  non¬ 
linear  vibratiems  of  sandwich  plates  are  discussed. 

Thus,  by  tbe  use  of  tbe  procedure  just  outlined,  we 
are  misbled  to  derive  tbe  appremimate  solutiem  to  an 
elMticity  problem  from  a  unified  point  of  view  and  in  an 
iiitagrated  manner,  solely  on  tbe  basis  of  tbs  variational 


of  Motion,  nlthoo^  aneik  tiMt  in  involved  nay 
lunre  been  wnll-bnown  na  iaolated  individaal  procedarea. 

the  rneiiniag  part  of  thin  paper,  ttie  nonlinear  vi- 
braticMia  of  hoRKyjeneona  »i»d  aandwicb  cylindrical  ahella 
will  be  inveati^ated  by  awaaa  of  the  propoaed  j^rocedare, 
with  a  ayaten  of  nonlinear  eqaationa  of  notion  of  cylin¬ 
drical  ahella  derived  in  the  firat  atep.  The  effect  of 
thidcneaa-ahear  defbmation  ia  inclnded.  The  problen  nay 
be  oonaidered  aa  an  extenaion  of  the  previoua  one  of  non¬ 
linear  vibration  of  ^t^ndwi^  plat>'u  [9],  and  the  reavlta 
arc  alao  re<hicible  to  eone  of  the  linear  reaulta  obtained 
previonaly  for  hoatogeneona  [12]  and  aandwich  cylindrical 
ahella  tdl. 
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Wpnlin— r  Boo>tlon»  of 
^fa«oq»n«om  and  Sandwich  CyltnJrieal  Shells 

BqiMtions  of  tli*  saadwidt  cylladrieal  ahall  wim>e 
dariirad  firat.  thoaa  of  tha  hoacganaona  ahall  will  than 
ha  ohtaiaaitla  aa  a  apaeial  Halting  caaa  by  patting  aqaal 
to  aaxo  tha  thidcaaaa  of  tha  faca  layara  of  tha  aandwidt 
ahall.  lha  cylindrical  ceordinataa  x,  a  •  aS  and  r  ara 
ehoaan  to  ha  in  tha  longitudinal,  cixcunfarantial.  and 
radial  diraetlotai.  raapactiwaly.  of  tha  ahall,  whoaa 
alddla  aorfaca  haa  tha  radiaa  a.  Tha  niddla  anrfaca  ia 
farther  daaignatad  aa  c  «  o  ao  that  tha  ralation  hatwaan 
tha  two  varitiblaa  r  and  a  ia  givan  by 

r  a  a  4-  a 

In  tha  a^iraction  tha  thicVnaaaaa  of  tha  inner  faca.  cora. 
and  outer  faca  layara  of  tha  ahall  axtand  froci  •h  to  -h^, 
•hi^  to  h^,  and  h^  to  h.  reopactively.  m  the  x-diraction 
tha  ahall  axtanda  froai  xaotox«l.  In  tlw  s-uiraction 
tha  ahall  ia  a  cloaed  ona  ia  tha  apacific  vibration  prob« 
lava  to  ha  diacaaaed,  but  tha  aquationa  darived  will  alao 
ba  applicabla  to  opan  ahalla.  In  particular,  tha  appropri¬ 
ate  boundary  oonditiona  for  an  opan  edge  aaeonatwit  aa 
wall  aa  thoaa  for  x«oonatant  will  ba  fomulatad. 

In  tha  caaa  of  mwII  dafOnaationa  and  anai-l  anglaa 
of  rotation,  aawll  in  a  aanaa  nodi  aa  apacifiad  in 
Movoahilov'a  book  [13],  tha  variational  aquation  of  notion 
ia  tha  nonlinear  theory  of  alaaticity  nay  ba  written  in 
cylindrical  ooordinataa  for  tha  aandwich  ahall  aa  followat 


wiwr*  and  ar«  two  arbitrary  aomanta  of  cina  and 
I  ^  T-^ii  *^r;  *  Txr.‘  ^s; 

=»  r^s.  ^  ^*V;  ~ 

-  rrri  ♦  ta;;  ~  <^«'  t*^si 

♦“ 

<T^*  a  (^*  ♦  Cjii;  «^ri  •“ 

,1 

'  ^ 

X^r,'  »  tgi  ^  T*-  -  tji^;  i»^,' 

^»X»'  *  XTra**  ■*“  Xrj,»  4^/ 

»  T«r  ♦  T„'  »>„•  -  tri  «j»- 

4f  *  ♦  fni  ‘*>>i  *  W«^/ 
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In  addition,  the  subscript  i  =  1,  2,  or  3  refers  to  the 
core,  inner  or  outer  face  layer  of  the  swdwich, 

Tjjsi,  are  the  stresses,  uiri  the  angles 

of  rotation  id)out  the  x,  s,  r  directions,  {^i,  Pri 

the  prescribed  surface  tractions  in  these  directions,  v 
refers  to  the  external  normal  direction,  ud  P;  is  the 
density.  While  the  volume  integral  in  the  equation  is 
to  cover  the  volume  of  each  of  the  three  layers,  the 
surface  integral  will  cover  only  those  portions  of  the 
surfaces  of  the  layer  on  %diich  tractions,  but  not  dis¬ 
placements.  are  prescribed.  The  equations  of  motion  and 
the  appropriate  boundary  conditions  of  the  sandwich  cy¬ 
lindrical  shell  will  be  derived  from  Eq.  (1)  by  carrying 
out  explicitly  the  integration  "ith  respect  to  z.  we  note 
here  that  this  will  involve  exactly  the  same  manipulation 
as  in  the  corresponding  linear  case  if  the  fogies  of  ro¬ 
tation  ate  independent  of  z,  %»hich  as  we  shall,  see  %fill 
be  assumed. 

As  in  reference  6  the  flexural  rigidit's-fS  of  the 
face  layers  will  be  neglected,  and  the  displ.’.ceiuents  are 
assumed  in  the  form 

U,  S  ut  i  Y  ,  Wj  s  ti  f  h,  V 

V/,r  Vf  Vj  *  V  >  (p 


The  angles  of  rotation  Wjjj  and  in  the  ccic  are  then 

a>  -sfSt  ^  ^  -  Yj 


Since  the  face  layers  have  been  taken  as  membranes,  their 
angles  of  rotation  about  the  x-  and  s-directions  are 
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The  angles  of  rotation  about  an  axis  normal  to  the 
shell  middle  surface  are  in  general  much  smaller  than 
i»jti  and  «»si  about  axes  lying  -*n  the  surface  and  are 
usumed  to  be  negligible.  Among  all  the  angles  of  ro¬ 
tation,  therefore,  only  y»xx  depends  on  s.  To  simplify 
the  formulation  of  the  problem,  this  z-dependence  is 
suppressed  and  is  assumed  to  take  the  simplified 
form 


since  the  original  z-dependence  of  at^i  de- 
cre2ise8  with  the  thic.  .iess-to-radius  ratio  of  the  shell 
and  lUxi  effects  only  the  nonlinear  terms,  the  simplifying 
assumption  should  not  introduce  much  inaccuracy  no  long 
as  the  shell  is  thin  and  the  nonlinearity  s<nall.  It  is 
also  noted  that  the  assumption  will  not  af  ftiffeet 
symnetrical  vibrations  which  are  to  be  disc-  ssed  later 
in  this  paper.  To  summarize,  the  angles  ot  rotation  now 
take  the  form 


Xj  ^  j  I 


w.,  -  u>i,  = 


- 


(4) 


iaJ,.  =  tc>ri  ~  uJf'j  0 


when  Eqs.  (3)  are  substituted  in  Eg.  (1)  and  inte¬ 
gration  is  carried  out  with  respect  to  z,  the  volume 
integrals  in  Eg.  (1)  yield  a  double  integral,  involving 
the  values  of  the  stresses  ^rsi,  "^rl* 

curved  boundaries  of  the  shell.  The  surface  integrals  in 
Eg.  (1)  yield  three  parts  after  integration.  The  first 
part  is  the  result  of  integration  over  the  curved  st'.r faces 
of  the  shell,  which  is  also  in  the  form  of  a  double 
integral,  when  only  tractions  are  prescribed  on  the 
curved  surfaces,  the  two  double  integrals  deduced  from 
the  volume  and  the  surface  integrals  may  be  condiined,  and^ 
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by  equating  to  zero  the  coefficients  of  6u,  6v,  6w,  6y,  ficD 
in  the  resulting  integrand,  the  following  stress  equations 


of  motion  of  the  shell  are  obtained: 


^ ,  ij’l  ‘'i)  * 


«)S 


(5) 


*(.!'>) ►V- ^3 

+  )v  h,(  1+ j)  -  ii-(eg*2-fj<,\)f=0 

-a4)J  t  (/.'  -  V*J]  -  ) 

K^-l I*b-  'h»'l>-i)- p, f 


where  ^J/c  A;,,f 

Eqs,  (5)  contain  '^nly  the  surface  tractions  pj*  ••• 
but  no  longer  the  boundary  values  of  the  stresses  ''t'xi' 
“'^rsi*  ^ri*  double  integrals  deduced  from 

Eg.  (1)  are  left  uncombined,  the  integrand  of  one  will 
yield  shell  equations  which  still  contain  the  latter 
stress  values,  and  that  of  the  other  will  yield  the 
stress  or  displacement  boundary  conditions  fet  the 
curved  surfaces  of  the  shell. 

The  remaining  two  parts  derived  from  the  surface 
integrals  of  Eq.  (1)  are  results  of  integration  over 
sections  across  the  thickness  of  the  shell,  and  they  are 
in  the  form  of  line  integrals  along  the  edges  x«constant 
and  s«constant.  From  these  line  integrals,  the  following 
appropriate  boundary  conditions  are  obtained* 
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Along  edge  x  =  constant, 

N^=  (^p.l  K-i)  I** 


<r- 

♦  ^ l< '’■J  Vi '-^ » 

-h.f"^'f^li-‘iA‘  "''V’V  '‘’ 

Mi,  ^  '  i:;  >s  >  ( •■'  i  >'*  ’  *■'"£  ^  ‘ 

Aong  an  edge  s  =  constant  (for  an  open  ^hell) . 

J-*, 

.  Cy  ^ 

I  n>  -h 

-h.  >^,  '• 


cv  \y  -  \^ 


jjy  vV= 


u=  u 
or  'Jr'^ 


■Jr  lV  ■ 


In  Eos,  (5)  to  (7)  overbars  in  general  denote  pre¬ 
scribed  quantitites.  Thus  ^  ,  ^s'  ^  **■* 

scribed  tractions  at  the  outer  boundary  *  -  'h.  Pj, 

those  at  the  inner  boundary  *  -  -h  and  Px.Ps,  Pr  ^^o»*  _ 
across  the  thicVness  of  the  shell.  biVevise.  u.  v.  w.  i'.  - 
are  the  prescribed  values  of  the  shell  displacements.  The 
shell-stresses  introduced  in  these  equations  are  defined 

by 
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(  m/,  , 

^  '  J 

(.Mil  j 

J 

( (TsT, 

(<r/ 

,  .  * 

)  ~  j 

( <^5.- 

)  = 

f  (  lJ,  . 

1 

1  )  - 

*.  1  )  d± 

r\  f  ^ 

where  the  inlegiations  are  understood  to  cover  ur.o  thick¬ 
ness  of  the  layar. 

The  stress-strain-displacement  relations  of  the 
shell  will  next  be  formulated,  tinder  the  assumptions  of 
small  uuyi  and  tuni  and  rero  "xi  strain -displacement  ’-c 
latinns  in  the  nonlinear  theory  of  elast>..  ity  are  [13] 


^  ^  J 
ax  ^  " 


63;  .  ^  ^  ^  ^ 


(  .  -  ^  diir  r 


>r. 

c 

-St^i 


* 

i  VX 

im 

r 


%rtiich  yield,  by  virtue  of  Eqs.  (3), 


1- 
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V'  ■■  3> 

-  f  (fe  ^  '  II-  - 

^xr-,  -  Y  fS 
^sr.  =  ‘f-  ^—r'  *  r 
e,.,  6,,  -  fi  ^  ^  f{-  -  i  ‘-si 
^SV  «S.  =  T  (If  f  ^  *f;  *  T  ‘  t 

^  ^  /  ^  qp  i,  ^  ?  I,  -  uJy, 

V>i  7L7s^*’'T5j  ^'' 


(9) 


where  ivi.  ‘J^’sx.  'yx2.  ''‘s2  are  given  in  Eg.^.  (4) 

The  stress-strain  r-'*  » Lions  for  the  laterials  of 
the  various  layers  arc 


iT,.  = 

(  <>/  r  K‘  €f,: ) 

^  - 
^  • 

t 

^xv  " 

M. 

^,,1  •  '■•y 

^$r/  -  ,«i( 

where  *  ^>/ ( i ' ‘^i  V-  A.=  /<?(  i t »>',•_)  .  e^  is  Young’s 

modulus,  Poisson's  ratio,  and an  additional  shear 
modulus  of  the  orthotropic  core.  Successive  substitution 
of  Eqs.  (9)  in  Eqs.  (10),  (2),  and  (R)  yields 
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A/s,'  - 

^  ‘^x: 


a/,f 

Aisxi  "  A.1sxi  -(3s, ^>1 


0.^;  =  <3^;  *3jrji  -  <3.5, 


where 


w..  --j  r..0^i ;  <.,/ « r  (If 

S  =  Jfs/  ■(  i  •  ?  2,  li,[  Ul,:  .  ,v,,-j 

N.  i  =jV.,(l%-j‘«<  --  ?A..,[  |i  ,  ..  a,;,,  ‘U„j 

N!,|=  fS:,i  ‘>i  ■  '0„  -Jill 


J  ^'1 


'^Xl  -  f  <•  XU  (  t  '  lx,  xi^i  It,  (  I,''  * 

Cj.;,  -  J  T,„  d  t  •  >  X,A„ll,  (‘f*  X) 

Ms„=Jrix,  ^  ^  -  a  Itj 

+^)l 


If 


As  those  1.'!  Eqr .  (9).  the  liwito  of  the  integrations  in 
Eqs.(12)  cover  the  thicVi-t  ss  of  the  liyer.  In  carrying 
out  these  integrations  the  tn:  n-  nhr  ll  jinpLio  of 

I  h2S  been  made  use  of  wherev  -r  applicable. 

In  the  expressions  of  a  shear  coclf  i  cient  Xj 

has  been  inserted  in  the  f.c..-.c  manner  as  ir  the  linear  case 
ri4,  15,  anv  ray  also  rc  de! ermir.nd  sir.rlarly.  The 
transverse  shear  forces  0x2,  0s2.  0x3,  ''’'•3  tero  be¬ 

cause  the  face  layers  have  been  assumed  to  be  membranes. 

Substiti’tion  of  Ees.  (11),  and  (12)  in  Cgs.  (5)  leads 
immediately  to  the  displacement  equations  of  motion,  al¬ 
though  the  results  will  not  be  recorded  heic. 

The  system  of  Eqs.  (5),  (€) ,  (7),  (11),  and  (12)  may 
be  readily  reduced  to  previous  results  fer  .'impler  special 
cases.  Thus,  by  letting  the  radius  a  equal  to  infinity, 
the  nonlinear  equations  of  sand'vich  plates  of  reference  9 
are  obtained,  cn  the  other  laud,  if  the  nonlinear  effect 
is  suppressed,  \/e  arrive  -:,t  one  of  the  simpler  systems  of 
linear  eoua^ions  of  sandwich  cylindrical  shells  of  reference 
P,  which  are  further  icducible  to  the  equations  of  homo¬ 
geneous  shells  of  reference  12. 
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Modification  auid  Simplification  of  Ecmations 


As  in  reference  5,  a  second  system  of  equations  of 
somewhat  better  accuracy  may  be  derived  by  assuming,  in¬ 
stead  of  Eqs.  (3), 


U,  =  ^  ^ 

V,  ^  \yf  if 


Uj  -  M.f  ;  f 

Vi  ■■  V'  7  )  f/' 


W,  -  iVj.  -  lA^  -  Ia/ 


The  results  are  similar  to  EqSi  (5);  (6',  (7^,  (11).  and 
(12).  In  fact,  they  are  obtainable  from  these  equations 
by  replacing  hj^  by  hj^  +  ho/2  in  the  last  two  of  Eqs.  (5) 
(except  the  h^  in  2hi3/3,  in  the  last  two  of  Eqs.  (6)  and 
those  of  Eqs.  (7)  (except  the  hj^  in  the  Units  of  the  inte¬ 
gration)  ,  and  in  t’  ■  expressions  of  Sx2.  *’,c3»  ^82.  ”s3' 
Nxa2*  ^v.s3<  '"’sx2-  Vsx3- in  mentations  (12).  vrith  the 

modifications  ..lade,  we  shall  designate  the  newly  obtaine-’ 
second  system  of  equations  as 

Eqs.  (5).  (6).  (7),  (li),  (12)  modified 

(5‘).  (6*),  (7’).  (11*). ^12') 

As  may  be  seen  from  Eqs.  (II)  and  (12),  by  having 
ta)<en  into  consideration  the  rotations  and  'Ogx 
core,  not  only  the  membrane  forces  H^l,  *^sl,  *^X8'',  **'** 

Nsxi  themselves  arc  effected,  but  also  the^  -ire  augmented 
by  the  transverse  shear  forces  Qxi  and  Conversely, 

the  transverse  shear  forces  are  also  augmented  by  the  mem¬ 
brane  forces.  Since  for  relatively  low  frequencies  the 
motion  of  the  shell  is  predominemtly  transverse  in  nature, 
the  contribution  of  Q^i  and  0,j  to  the  membrane  forces 
should  be  of  less  importance  than  the  contribution  of  Nxi, 
**sl,  imd  N,j{i  to  the  transverse  shear  forces.  The 

contribution  of  and  Qg^to  the  moments  is  also  small, 
for,  according  to  Eqs.  (ll),  they  are  multiplied  by  both 
a  small  angle  of  rotation  («xl  or  u)*].  and  a  small  factor 
hj^^/3a.  As  a  simplification,  the  contributions  of  Q^i 
and  Qgj  to  the  membrane  forces  and  moments  are  neglected 


in 

in  Eqs.  (11).  If  we  are  primarily  interested  in  low  fre¬ 
quencies  for  %ihich  the  motion  xs  predraninently  transverse, 
the  inertia  terms  involving  u,  v,  y  emd  v  become  of  much 
less  importance  than  the  tramsverse  inertia  term  involving 
w  and  may  *lso  be  neglected.  Incorporating  both  types  of 
simplifications,  we  shall  designate  the  results  obtained 
from  the  first  system  of  equations  by 

Eqs.  (5),  (6).  (7>,  (11).  (12)  simplified 

(Sa),  (6a).  (7a).  (11a),  (12a) 

The  second  system  of  equations  may  similarly  be  simplified 
and  denoted  by 

Eqs.  (5*).  (6‘).  (7*:  .  (U*).  (12*)  simplified 

(5a').  (6a'),  '.7?.').  (11a').  (12a') 

For  sandwich  shells  with  soft  cores  membrane  forces 
and  moments  in  the  core  may  further  be  ncg’,?cted,  that  is, 
may  put 


/\i^i  =  A/', I  ■“  A'x./  ■  U;  ~  —  0  (13) 

Contrary  to  this,  we  have  in  the  case  of  hoi!K>geneous  shells, 
for  which  h2  ■  0, 


=  'Vsr  =  O  (14) 


In  the  axially  syrmietrical  case  with  sero  surface  tractions, 
Eqs .  ( 5a ' )  become 


where 
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-  2  e,h,  I  u'*  ^  y ^ 

<\)5i  =  +  v/^v  ^ 


■  ’  !>i  -  '^52./ 


with  primes  denoting  differentiation  with  respect  to  x. 

The  final  simplifications  given  by  Egs,  ;1  ')  aoid  (14)  will 
be  it-:  irporated  in  Egs.  (15)  in  the  next  S'rctions  in  which 
axially  syTtcne*  rical  vibrations  of  sandwicn  and  homogeneo’'r 
cylindrical  shells  are  discussed. 


nonlinear  Vibrations  ot  Horoogeneous  and 
SdPdvdrh  0>lin»^r~cal  Pbc-lls 
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Tt  should  be  emphasized  at  this  point  that  each  of 
the  above  systems  of  shell  eqi'ations  is  essentially  the 
result  of  putting  to  zero  the  coefficients  of  the  varia¬ 
tions  of  the  displacements  in  the  integremds  of  the  now 
simplified  veuriational  equation  of  motion,  simplified  in 
that  it  no  longer  contains  the  integration  with  respect 
to  z.  For  those  of  the  shell  equations  that  cw  be  solved 
2md  satisfied  exactly  in  a  given  problem,  the  correspond¬ 
ing  --.'.efficients  will  simply  drop  out  of  the  simplified 
variational  equ  -.tion  of  motion.  On  the  oth*»r  hand,  for 
those  shell  equations  that  cannot  be  solved  exactly,  the 
corresponding  coefficients  will  remain.  Subsequent  vari¬ 
ational  aj^roximation  may  then  be  performer,  %diich  maVes 
the  latter  ecoiations  eventually  also  satisfied,  at  least 
approximately  in  the  variational  sense.  This  procedure 
will  now  be  demonstrated  by  the  following  .itscussion 
axially  symmetrical  vibration  r*  homogenecus  and  sandwich 
cylindrical  shells  with  itmovable  hinged  edqes. 

The  homogeneous  cylindrical  shell  will  be  considered 
first,  for  which  Eqs.  (14)  to  (16)  yield 


i  y  (  f"'*’  i  2.X^A  ( Q 


(17) 


(18) 


where  the  only  rotation  component  involved  is  now 


'  i  (  '^) 
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The  subscript  1  has  been  dropped  from  all  notations  in 
Eqs,  (17)  to  (19).  The  shear  modulus  u  is  now  for  am 
istropic  material,  and  the  shear  coefficient  has  the 
usual  value  of  n2/12  ri4l,  although  it  may  also  be  as¬ 
signed  the  value  of  infinity  for  the  purpose  of  suppress¬ 
ing  the  transverse  shear  effect. 

From  Eqs.  (17)  and  (10)  u"  is  first  eliminated.  Into 
the  result  are  s^ibstituted 


^v=  W  Tit)  ^  cos  ^  (20) 


which  satisfy  the  conditions  of  sero  deflection  amd 
moment  at  the  hinged  edges  x  =  O.JI;.  It  xs  then  found 


that 

5 


_ I  _ _ 


(21) 


In  Eqs.  (20)  .ud  (21)  we  have  > =  nnh/i  ,  with  n  *  1,2,3,  , 
designating  the  number  of  naif-waves  in  the  l»n'’th  ^  or  the 
shell,  and  ■'■(t)  is  the  unknown  time  function,  u  is  next 
solved  from  Eq.  (17).  into  \d\ich  w  and  Vai?  substituted. 
Together  with  the  use  of  the  boundary  conditions  of  u  =»  0 
at  the  iiranovable  edges  x  *  0,  ji  ,  we  find 


Eqs.  (20)  to  (22)  thus  satisfy  exactly  all  the  boundary 
conditions  and  the  governing  Eqs,  (17)  and  (IP).  However, 
it  is  easily  verified  that  they  do  not  satisfy  the  remain¬ 
ing  governing  Eq.  (15).  Since  the  left  side  of  Eq.  (19)  is 
the  coefficient  of  6W  in  the  double  integral  (actually  re¬ 
duced  to  a  single  integral  of  x  in  the  present  plane-strain 
problem)  in  the  variational  equation  of  motion,  and  since 
we  now  have 


r 


K> 


=  S-i  ^  6(  V\J-c) 


2^ 

according  to  Eqs.  (20)  we  may  carry  out  explicitly  the 
integration  with  respect  to  x  over  the  length  of  the 
shell  amd  put  the  coefficient  of  6(Wt)  equal  to  zero. 

Thus,  there  results 


■c  X  t  >‘vt:  oii  r  ^  =  O 


(23) 


This  in  essence  has  completed  a  second  vari  *tior<al  approxi¬ 
mation,  with  respect  to  the  second  amd  last  remaininq 
space  coordinate  x,  although  the  approx irriot< on  applies 
to  only  one  of  the  governing  differential  e-^jations.  If 
needed,  the  procedure  could  have  also  been  applied  to  any 
other  or  all  of  the  differential  equations  iwd  lx>undary 
conditions. 

To  determine  the  nonlinear  freqxjency  it  is  convenient 
to  carry  out  a  last  variational  approximation  with  respect 
to  time,  in  connection  with  Eg.  (23) .  The  left  side  of  the 
eouation  is  essentially  the  coefficient  of  6(wt)  in  the 
variational  equation  of  motion,  which  has  now  been  much 
simplified,  since  the  time  integral  is  the  only  one  re¬ 
maining.  The  integration  with  respect  to  time  may  be 
carried  out  explicitly  by  first  assuming  say,  t  »  slnu't, 
and  by  selecting  tg  ■=  0  and  tj^  »  2TT/tt)  as  the  limits  of 
integration.  The  coefficient  of  4w  is  finally  put  equal 


to  zero  to  yield 
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to 


It  is  interesting  to  note  that  the  approximate  nonlinear 
frequency  is  independent  of  02* 

In  a  similar  manner  we  may  investigate  the  vibration 
of  the  sandwich  cylindrical  shell,  to  which  Egs.  (15)  and 
(16)  together  with  the  simplifications  in  Eqs.  (13)  are 
applicable.  The  equations  of  motion  in  tetms  of  the  dis¬ 
placements  are  of  the  final  form 

+  r  ^ 

tv  iv  tc 

vv".! 

,  V  iS) 

2  )ij.|  a't  1/  y '  t-1  kV'*"  [  w  " 

”  ^  I  ^  It'  -<•  -I- w'  ‘j  -  c) 


The  shear  coefficient  may  be  taVen  eoual  ‘o  I  for 
ordinary  sandwich  shells  fls,  h],  although  ic  may  also 
be  set  equal  to  infinity  for  the  purpose  of  suppressing 
the  transverse  shear  effect. 

Eqs.  (25)  are  entirely  similar  to  Eqs.  (17)  to  (19) 
and  the  seune  method  of  solution  is  applicable.  The  results 
are 


w  r  t  S'*' 

Y  -  -  vAy  r  ,  Y  = 
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The  nonlinear  freouency  is  again  given  by  {he  approximate 


expression 


which  is  now 


i-y  -  5^,  +  2?  ^3 


Oj'"- 


SJZJl*! _ L 


r,>yx, 


2  4^^ 


(26) 


Eg.  (26)  is  reducible  to  the  result  for  th-:  nonlinear 
frequency  of  sandwich  plates  (9]  by  letting  a  equal  to 
infinity,  and  to  the  result  for  relatively  low  linear 
frequencies  of  sandwich  cylindrical  shells  (81,  by 
dropping  the  nonlinear  term  (3/4)  (v?/2h)^. 

The  frequencies  of  nonlinear  symmetrical  vibrations 
of  homogeneous  and  sandwich  cylindrical  shells  have  thus 
been  determined  solely  on  the  basis  of  the  variation  e- 
quation  of  motion,  by  carrying  out  variational  approxi¬ 
mations  wherever  needed.  Further  refinements  may  be  made 
in  the  analysis  and  the  discussion  may  readily  be  extended 
to  non-symmetrical  vibrations  and  to  other  types  of 
boundary  conditions,  but  the  relatively  simple  cases  pre¬ 
sented  here  are  clearly  sufficient  for  illustrating  the 
variational  procedure  that  is  being  proposed. 


Discussion  of  Results 


The  results  in  Eos.  (24)  amd  (26)  give  the  freouencies 
for  the  lowest  faunily  of  axially  synuaetrical  modes  of  vi¬ 
bration  of  the  cylindrical  shells,  which  are  predominate¬ 
ly  transverse  in  nature.  The  effect  of  thiclcness  shear 

deformation  is  associated  with  the  ^  -  or^'j^-  terms,  and 

2 

its  importance  is  seen  to  vary  directly  with  A  ,  Since 
putting  X.  orX'i  equal  to  infinity  is  equivalent  to  the 
suppression  of  the  shear  effect,  Eos.  (24)  and  (26)  show 
that,  disregarding  nonlinear  aind  curvature  effects  at 
this  moment,  the  shear  effect  by  itself  becomes  negligible 
if 

^  ~  ^ )  (homogeneous) 

or 

>  ~~ —  (sandwich) 

’t  > 

These  conditions  apply  to  plates  as  wll  ns  shells.  Since 
r2ri|  for  ordinary'  sandwic’.  structures  is  <  sually  of 
order  of  between  10  amd  100.  c»»e  shear  eitect  is  much  more 
important  for  sandwich  than  for  homogeneoi.e  plates  and 
shells.  In  general,  the  shear  effect  should  be  considered 
for  sandwich  structures  [16]. 

The  and  h^/a^-  terms  in  Zqs.  (24)  and  (26) 

reflect,  respectively,  the  nonlinear  and  curvature  effects, 
which  are  completely  uncoupled  from  each  r. Uier,  Tents 
associated  with  the  coupling  between  the  two  effects  would 
be  present  in  the  equations,  if  we  had  employed  more  exact 
expressions  of  the  rotation  components.  Since  it  is  the 
approximate  expt-essions  of  the  rotation  components  in  Eqs. 
(4)  which  have  been  used  and  vd’.ich  are  the  seune  as  those 
for  plates,  the  nonlinear  terms  in  Eqs.  (2^*)  and  (26)  are 
also  essentially  the  same  as  those  for  plates.  In  the  case 
of  sandwich  plates,  the  nonlinear  effect  was  discussed  be¬ 
fore  in  reference  9,  where  it  was  found  that  the  nonlinear 
effect  would  overshadow  the  shsar  effect  if 


(27) 
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